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ABSTRACT 

Taking three independent approaches, we investigate the simultaneous con- 
straints set on the cosmic star formation history from various observations, in- 
cluding stellar mass density and extragalactic background light (EBL). We com- 
pare results based on: 1) direct observations of past light-cone, 2) a model using 
local fossil evidence constrained by SDSS observations at 2 ~ (the 'Fossil' 
model), and 3) theoretical ab initio models from three calculations of cosmic star 
formation history: (a) new (1024)'^ Total Variation Diminishing (TVD) cosmolog- 
ical hydrodynamic simulation, (b) analytic expression of Hernquist & Springel 
based on cosmological Smoothed Particle Hydrodynamics (SPH) simulations, 
and (c) semi-analytic model of Cole et al. We find good agreement among the 
three independent approaches up to the order of observational errors, except 
that all the models predict bolometric EBL of /tot — 37 — 52 nW m~^ sr~^, 
which is at the lower edge of the the observational estimate by Hauser & Dwek 
(2001). We emphasize that the Fossil model that consists of two components — 
spheroids and disks — , when normalized to the local observations, provides a 
surprisingly simple but accurate description of the cosmic star formation history 
and other observable quantities. Our analysis suggests that the consensus global 
parameters at z = are: = 0.0023 ± 0.0004, Jebl = 43 ± TnWm-^sr'^ 
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1. Introduction 

Targets of cosmological investigations during the last decade have largely shifted from 
the global, geometrical properties of the universe to the detailed contents thereof, particularly 
the origin and evolution of structure from smooth initial conditions. Of the components 
added since the decoupling of matter and radiation, stars are the most easily observable 
feature (emitting the most electromagnetic energy) and of the greatest historical significance. 
Approximately 6 ± 2% (e.g., Fukugita & Peebles 2004) of the baryons of the universe have 
condensed out into stars at the current epoch, and the time history of the process whereby 
this occurred is the subject of intense study at present. What is not widely realized is that 
there are three independent approaches by which this study can be pursued: 

1. We can use the universe directly as a time- machine: in areas of the sky like the Hubble 
Deep Field (Williams et al. 1996), attempts have been made to analyze the past light- 
cone. The rate of transformation of gas into stars is measured as a function of the 
look-back time or redshift, which produces what has been called the "Madau Diagram" 
(Madau et al. 1996). There are three major uncertainties involved in this diagram: (1) 
dust obscuration; (2) relation between emission and star formation rate (SFR), which 
includes the uncertainty in the initial mass function (IMF); (3) faint-end slope of the 
luminosity function, especially at high- redshift. All of these three contribute roughly 
on the same order. Recent observational estimates at 2; > 4 suggest a decline in SFR 
with increasing redshift (Iwata et al. 2003; Bunker et al. 2004; Bouwens et al. 2004; 
Giavalisco et al. 2004; Ouchi et al. 2004a; Bouwens et al. 2005), but due to uncertainties 
in dust extinction for ultra-violet (UV) light and the faint-end slope of the luminosity 
function, the trend does not seem to be inconsistent with being constant at 2; > 3 (See 
also Figure 2 of Nagamine et al. 2004). 

2. The second method predated the direct method outlined above. It uses the fossil 
record from our own and other galaxies to catalog stars of different ages and to, by 
an essentially demographic investigation, unravel the history of star formation in the 
local Universe using the theory of stellar evolution. When direct age estimates are not 
available, as is the case for external galaxies, color distributions can provide a useful 
proxy to their history, distinguishing between systems with and without ongoing star 
formation activity, although there are degeneracies among age, amount of stars formed, 
metallicity, and IMF. This method was pioneered by Searle, Sargent, & Bagnuolo 
(1973). Here one makes a simple ansatz for the star formation history — chosen in 
the quoted paper to be as a simple declining exponential — and an IMF, and then 
computes, via the standard theories of stellar evolution and stellar atmospheres, the 
evolving spectral output from the assemblage of stars. They found that one could easily 
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and plausibly simulate the observed colors of different types of stellar systems. We will 
update the quantitative input to this theory, leaving the essential basis unchanged, 
and find that this "Fossil" model fits with uncanny accuracy (given the simplicity of 
the method) modern observations as obtained by the Sloan Digital Sky Survey (SDSS; 
York et al. 2000). We will then use SDSS results to determine the parameters of the 
Fossil model. 

3. The third method is to attempt to solve the problem by ab initio theory One takes the 
initial conditions from the observed cosmic microwave background radiation, adopts a 
favorite cosmological model, e.g., the standard concordance A cold dark matter (CDM) 
model (Ostrikcr & Steinhardt 1995; Perlmuttcr ct al. 1998; Ricss ct al. 1998; Spcrgcl 
et al. 2003; Tcgmark et al. 2004), combines this with the standard physical equations 
of gravity, hydrodynamics, atomic physics, radiative transfer, etc., and simulates the 
history of the universe by numerical integration, obtaining as the most important by- 
product the history of star formation (e.g., Nagamine et al. 2000, 2001b; Ascasibar et al. 
2002; Springel & Hernquist 2003b; Nagamine et al. 2004; Nagamine et al. 2006). One 
could also employ the so-called semi-analytic models of galaxy formation (e.g., Kauff- 
mann et al. 1993; Cole et al. 2000; Somerville et al. 2001), where a dark-matter halo 
merger-tree, cither constructed using the extended Press-Schechter formalism (Bond 
et al. 1991; Lacey & Cole 1993) or N-body simulations, is supplemented with a set of 
equations to model transformation of gas into stars within virialized dark matter halos. 

If the observational errors are too significant or the theoretical modeling in error, these 
three independent approaches might not agree with one another when they should. We need 
to pay particular attention to the errors that are involved in various observations. The quan- 
tities at near zero redshift are usually measured without too serious systematic errors. On 
the other hand, high redshift observations, such as those of star formation rates are subject 
to large systematic uncertainties from the obscuration and the faint-end slope of luminos- 
ity functions. Therefore, we take the attitude that we constrain the three approaches with 
zero-rcdshift quantities, while high-redshift observations are used to examine the consistency 
within allowed errors. We also note that there are some uncertainty such as IMF, but these 
uncertainties are known and we can see how they propagate to our calculations or even con- 
strain those uncertainties. With understanding those uncertainties, we can conclude that the 
agreement of the three approaches is in fact good. In addition, there is an independent check 
of all of these approaches. Each approach predicts two numbers: (1) the local bolometric 
electromagnetic energy output per unit volume jhoii^ = 0) [erg sec~^ cm"''] and (2) the local 
energy density U{z = 0) [erg cm~^], which is measured as the extragalactic background light 
(EBL). The first number jboi(^ = 0) provides a measure of the current star formation rate. 
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while U{z — 0) constrains the integral of all past star formation. The local energy out- 
put provides the normalization of the first two quoted approaches and consequently cannot 
be used as an independent check, but the EBL provides an independent test of the three 
approaches. 

The paper is organized as follows. In Section 2 we describe the fossil approach character- 
ized by delayed exponential decay models of cosmological star formation and investigate the 
resulting color and mean age as a function of characteristic decay time and metallicity. As a 
result, we construct the 'Fossil' model using two components, bulge and disk, based on the ob- 
served colors of galaxies. In Section 3 wc consider the star formation history, as predicted by 
the Fossil model and other ab initio theoretical models, and compare it with what is inferred 
from observation. In Section 4 we discuss the stellar mass density, a time integrated version 
of the star formation history, which is taken as a consistency test. We then compute the EBL 
in Section 6 to study whether what is expected from the star formation history is consistent 
with observations. Conclusions are given in Section 8. When necessary, we adopt a cos- 
mological model with parameters {QM,^A,^b, h,n,as) = (0.3,0.7,0.04,0.7,1.0,0.9), where 
h = -ffo/(100 km s~^ Mpc~^) is the Hubble parameter. Where h is not explicit, h = 0.7 is 
assumed. 



2. The Fossil model 

Scarlc. Sargent, & Bagnuolo (1973) showed that the colors of galaxies can be reproduced 
well by simple exponentially decaying star formation histories, p^, = A exp(— t/r), where r 
is the characteristic decay time. Here we revisit their model with a modern perspective. 
Throughout this paper, we take t = to be the epoch of the Big Bang. 

The original model of Searle et al. (1973) has the disadvantage that one has to assume 
the onset time of star formation by hand, increasing the number of free parameters. Therefore 
we adopt the delayed exponential model, 

= A(i/r)exp(-t/T). (1) 

This model has a desirable feature that the SFR vanishes at t = 0, and it has a peak at 
t = T followed by an exponential tail with a characteristic time-scale r. The model with a 
very small value of r resembles the instantaneous burst at t = 0, and the one with t tn 
(Hubble time) resembles the model with a constant SFR. The delayed exponential model has 
the same minimal number of free parameters as the Searle et al. (1973) model, a time-scale 
T and a normalization A. Good agreement with other approaches, as we will show later, 
justifies the use of this functional form for the Fossil model. 
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Figure 1 shows the delayed exponential models of star formation history as functions of 
cosmic time (panel [a]) and redshift (panel [b]) for different values of r as described in the 
caption. The normahzation A is adjusted so that the integral of the SFR, 



where pcrit is the critical mass density and to is the current age of the universe. The thick 
blue solid line is a composite two population model with r = 1.5 Gyr and 4.5 Gyr which will 
be described below. 

Wc note that the integral of SFR (Equation [2]) is not equal to the stellar mass density 
observed today because of gas recycling, i.e., gas lost by aging stars. The recycling fraction 
R, defined to be the ratio of the amount of recycled gas to the total amount of gas that was 
initially converted into stars, depends on the IMF and is a function of time. We show the 
recycling fraction in Figure 2 computed using the population synthesis model of Bruzual & 
Chariot (2003, hereafter BClibOS). This shows R ~ 0.45 at t = 13.5 Gyr for the Chabrier 
(2003) IMF integrated from 0.01 to 100 Mq^. In other words, approximately 55% of the initial 
gas remains in stars and stellar remnants (45% in shining stars and 10% in remnants, i.e., 
white dwarfs, neutron stars and black holes, for which BClib03 assumes the mass, 0.55, 1.4, 
and 2.0 Mq, respectively.) Therefore the normalization given by Equation (2) corresponds 
to the stellar mass density ili,{z = 0) 0.0022. We note that the recychng fraction of the 
SSP for the Salpeter IMF (0.1 - 100 Mq) is i? = 0.32 at t = 13.5 Gyr. 

We then compute the colors of galaxies for each model using BClib03 with the Chabrier 
IMF. The delayed exponential models give u — g and g — r colors as shown in Fig ure 3 a at 
to — 13.5 Gyr. The three lines stand for metallicities: Z/Zq = 0.2 {black open squares), 1.0 
{blue solid triangles), and 2.5 {red open triangles), with r varying along the line from top 
(0.1 Gyr) to bottom (10 Gyr) as indicated in the legend. Also shown as contours is the color 
distribution of SDSS galaxies taken from Figure 7 of Blanton et al. (2003b). 

The histogram (normalized arbitrarily) of the SDSS galaxy color projected onto each axis 
is shown in Figure 3b, with double Gaussian fits by eye. The two rectangular boxes delineated 
with dots indicate the locations of the Gaussian peaks with the width twice the dispersion: 
u — g = 0.94 ± 0.25 and g — r = 0.43 ± 0.14 for the blue peak, and u — g = 1.6 ± 0.25 and 
g—r = 0.88±0.07 for the red peak. Much has been said of the bimodal color distribution seen 
in the local galaxies from the SDSS (Strateva et al. 2001; Blanton et al. 2003b; Kauffmann 



^BClibOS integrates the Chabrier IMF from 0.1 to IOOMq. We applied the correction for brown dwarfs 
that reduces R from BClibOS by a factor of 0.94. 




(2) 
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et al. 2003; Baldry et al. 2004; Balogh et al. 2004; Brinchmann et al. 2004). The two distinct 
distributions correspond to non-star forming and star-forming galaxies, but roughly also to 
the redder spheroid-dominated component and the bluer disk galaxies. This is the concept 
first noticed by Baade (1944) in his proposal of the two population model. He also pointed 
out that the two types of stellar populations had been recognized as early as 1926 by Oort 
(1926). 

We represent the peaks of SDSS galaxy color distribution with two parameters (r, Z /Zq). 
Figure 3a shows that the models with r = 1.5 Gyr and 4.5 Gyr fit the colors of the spheroid 
and disk components reasonably well, respectively, and the variation of metallicity gives a 
further tuning. We find that it is difficult to realize the colors {u — g,g — r) = (1.6,0.9) 
with BClib03 even when we push the metallicity to Z/Zq = 2.5, which is unlikely to be a 
representative of spheroid component even if the most massive elliptical galaxies may have 
such a high metallicity. We present the stellar mass weighted mean metallicity for spheroids 
(early type galaxies and bulges of disk galaxies) of different morphological types in Table 4, 
as calculated in Appendix A. As a result, we construct a composite 'Fossil' model — we 
call it Fossil because it is based on the fossil evidence in the present universe, i.e., colors 
of galaxies — with the parameters (r [Gyr], Z/Zq) = (1.5, 1.5) for the spheroid component 
and (4.5, 0.8) for the disk component ^. These parameters are indicated by the asterisks in 
Figure 3b. The blue asterisk for the disk component is located at the center of the box, 
but the red asterisk for the spheroid is shghtly off to the corner. The mass- weighted mean 
metallicity of the spheroid component we obtain in Table 4 is 1.3 Zq, but in order to keep the 
red asterisk within the dotted box in Figure 3b, we choose to adopt 1.5 Zq for the spheroid 
component. 

Next, we determine the normalization of the two components. We choose to normalize 
each component by requiring a match of the energy output in K- and r-bands with observa- 
tions (see Appendix B). Our normalization corresponds to the bulge(spheroid)-to-disk stellar 
mass ratio (Mfc/M^) = 1.35, or equivalently Mb : Ma = 57% : 43%. An alternative method 
to determine the normalization uses the bulge-to-disk luminosity ratio and the mass-to-light 
ratio for each morphological type of galaxies, as outlined in the caption of Table 1; this gives 
(Mfo/Md) = 1.9, or M(, : = 66% : 34%. Given the uncertainties involved, we consider the 
bulge stellar mass fraction of 55 — 70% to be observationally allowed. For this work we adopt 
57% obtained in Appendix B, but the exact value is not important to the major conclusions. 



^Here we are representing the entire galaxy population with two delta functions on the color-color plane. 
It is possible in the future to extend Figure 3 into a three dimensional color space by adding another color 
as a z-axis, and assign a set of parameters (r, Z) to every grid cell in the 3-d color space. That would be a 
more smooth, continuous representation of stellar mass distribution in the color space. 
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Figure 4 shows the mass-weighted mean stellar age for the delayed exponential model, 
calculated as 



The mean age of disk stars (r = 4.5 Gyrs) turns out to be tage = 7 Gyrs, which is consistent 
with the observational estimate of 5 — 7 Gyrs (Rocha-Pinto et al. 2000; Robin et al. 2003; 
Naab & Ostriker 2006). We note that our estimate includes the stellar remnants. Had one 
omitted the contribution from stars that died in the history of the universe, then the mean 
age would be younger by approximately 0.8 Gyr, as indicated by the arrows. The main 
properties of our Fossil model at the present epoch are summarized in Tables 2 and 3. It 
may be helpful to note that the historical development of the present distribution of galaxies 
is of no import for the Fossil model. It makes no difference at all to the accounting if we had 
dry, wet, or no mergers. All that matters is that the local SDSS inventory of a representative 
volume enables us to estimate the age distribution of star formation. 



In Figure 5, we show the cosmic SFR density as a function of redshift. The blue solid line 
is the Fossil model described in the previous section, and the black short-dash long-dash line 
shows the result of a new Eulerian TVD hydrodynamic simulation. This simulation has a 
comoving box size of 85h'^ Mpc and 1024^ hydrodynamic mesh, and the code is similar to the 
one used by Cen, Nagamine, & Ostriker (2005) and Nagamine et al. (2004). The physical cell 
size is 83/(1+2;) kpc. The mean baryonic mass per cell for this simulation is nigas — 3.63 x 
10^ h~^MQ and the dark matter particle mass is uibm = 1-58 x IO^H'^Mq (512^ particles). 
The cosmological parameters are {Qu,^A,^b,h,n,(Tfi) = (0.31,0.69,0.048,0.69,0.97,0.89). 
The TVD result is extracted from the simulation without additional processing except for 
the boxcar smoothing over ±3 bins in the redshift axis. The amount of gas converted into 
stars is a direct output of the simulation, so there is no uncertainty as to stellar IMF and 
no freedom for the normalization. For the details on the star formation and supernova 
(SN) feedback prescription, we refer the readers to Cen & Ostriker (1993); Nagamine et al. 
(2001b); Cen et al. (2005). 

We also consider two other models: one from Hernquist & Springel (2003, H&S model) 
who derived an analytic approximation for the SFR history as a function of redshift based on 
their cosmological Smoothed Particle Hydrodynamics (SPH) simulations. The H&S model 
shown in green long-dashed line takes the form 




(3) 



3. Cosmic star formation history 



l + a(x-l)3exp(/3x^/4)' 




(4) 
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where x{z) = {H{z)/HqY/^. With star formation and feedback models adopted in Springel & 
Hernquist (2003a), the parameters are a = 0.012, f3 = 0.041, and p^(0) = 0.013 Mq yr'^ Mpc^^ 
The normalization is fixed to give the local SFR density (Hernquist & Springel 2003). 

The other model represented by red short-dashed line is the semi-analytic model GAL- 
FORM by Cole et al. (2000, hereafter SA model). The cosmological parameters that H&S 
and SA adopted (Qm, ^^a, ^b, h, n, ag) = (0.3, 0.7, 0.04, 0.7, 1.0, 0.9), differ slightly from those 
of the TVD simulation. We confirmed that this slight difference is not important for our 
analyses. Table 3 summarizes the basic characteristics of these models. 

For the observational data shown in Figure 5, we assume the Chabrier IMF (0.01 — 
IOOMq). The UV luminosity densities puv are converted into SFR by 

puv [erg Hz-i Mpc'^] = C [Mq yr'^ Mpc-^], (5) 

where the parameter C can be derived from BChb03 as 

C = 1.24 X 10^^ (6) 

for an exponentially decaying star formation history with r = 5 Gyr for the Chabrier IMF 
and solar metallicity. Here puv is computed by averaging the fiux over ±150 A centered at 
A = 1500 A at age t = 10 Gyr. The parameter C depends on the age of the stellar population 
and reaches a plateau after t ~ 100 Myr for r > 1 Gyr models. The value in equation (6) 
differs by a factor 1.6 from the value given by Madau, Pozzetti, & Dickinson (1998) who 
performed the same calculation for the Salpeter IMF (0.1 — 100 M©). 

The correction due to dust obscuration is important. For observations of nearby uni- 
verse, the Balmer emission lines (Gallego et al. 1995; Tresse & Maddox 1998; Hopkins et al. 
2000; Pascual et al. 2001; Tresse et al. 2002; Nakamura et al. 2004) can be used to estimate 
the SFR with reasonable rehabihty at z ~ 0. For higher redshift, where SFR is derived 
from the UV hght, estimates of extinction are more uncertain. Steidel et al. (1999) took 
E{B — V) — 0.15, somewhat smaller than the value inferred for the local galaxies. At some 
point at high redshift, the dust content should start decreasing as a function of increasing 
redshift, but whether this is seen towards higher redshift or not is controversial even today. 
Bouwens et al. (2005) suggest a modest decline in the dust extinction effect from z = 3 to 
z = 6, but Thompson et al. (2006) indicate the extinction correction of a factor of 5 at 2; = 6, 
which differs little from that ai z = 3 — A. In this situation we must be content with with 
our fiducial choice of the conventional prescription of Steidel et al. (1999) that the SFR are 
corrected by factors of 2.7 {z < 2) and 4.7 {z > 2), as also supported by the subsequent 
work by Reddy & Steidel (2004). We apply this correction to most of the data (after re- 
moving the corrections by individual authors) , keeping in mind that varying dust extinction 
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as a smooth function of redshift is needed to obtain more reliable estimates of SFR. Note 
that the faint-end slope of the luminosity function at high-redshift is poorly constrained and 
SFRs are derived under different assumptions (we do not dare to standardize them except 
for extinction corrections), which we review briefly in what follows. 

Steidel et al. (1999, filled pentagons at 2; = 3, 4) derived their UV luminosity density by 
integrating the luminosity function with the faint-end slope of a = —1.6 down to O.IL*. We 
take the SFR from their Fig. 9, and convert it to the A cosmology. Giavalisco et al. (2004, 
open triangles at 2; = 3 — 6) integrated the Schechter fit with a = —1.6 to 0.2L*. Ouchi et al. 
(2004a, filled circles) integrated the luminosity function with a = —2.2 to O.IL*. Bouwens 
et al. (2005, filled square at 2; = 6) integrated their Schechter fit with a = —1.74 to 0.04L*^3. 
Their sample goes to the faintest magnitude compared to others, yet the result still indicates 
some dechne in SFR density from 2; = 3t0 2; = 6bya factor of ~ 0.7. Thompson et al. 
(2006) performed the SED fitting to the 6 photometric broadband measurements by ACS 
and NICMOS of the Hubble Ultra Deep Field determining the extinction for each sources. 
They then utilized the star formation intensity distribution function (Lanzetta et al. 2002) to 
correct for the surface brightness dimming effect. Their result is consistent with the earlier 
result by Thompson (2003). Here we take the extinction uncorrected data listed in their 
Table 2 and applied our correction as mentioned above. Since their errors are large and the 
error bars overshoot to the outside of the plotted range of Figure 5 when put on our version 
of data points, here we omitted the error bars for this data. 

Extinction corrections are not applied to the data derived from X-ray (Norman et al. 
2004) and submillimeter observations Barger et al. (2000, taken from Table 2 of Hopkins 
(2004)), as they are less subject to dust extinction effects. The data by Heavens et al. (2004) 
at very low-redshift from the MOPED (Heavens et al. 2000) algorithm already include the 
correction for dust for each galaxy. The data of Nakamura et al. (2004) are derived from 
hydrogen B aimer lines including a consistent correction for extinction in individual galaxies. 
We note that uncertainties for z > 1 could be much larger than those represented here, while 
those for lower redshift are probably not too large since dust extinction is constrained well 
with the use of Balmer lines. 

With large uncertainties we discussed in mind. Figure 5 shows a good agreement among 
the three independent approaches: observed SFR, Fossil model, and ab initio models. All 
estimates agree with each other to within about a factor of two, which is also the size of the 
scatter in the data. We see some declining trend of the SFR in the SA model at 2; > 4, which 
contrasts to the other numerical ab initio models and the Fossil model that show a roughly 
constant rate of star formation. Our argument given here should be understood in the sense 
that if observations would receive too large systematic errors, the consistency with the two 
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other approaches is endangered. 

We also underline the conclusion from the Fossil model that the spheroidal component 
formed primarily at high-redshift {z > 1.5), and the majority of the disk stars formed at 
low-redshift {z <1). The sum of these two components falls in-between the curves of other 
models. 

A note on the uncertainty due to our poor understanding of star formation is appro- 
priate. The theoretical calculations are most secure in their predictions of the SFR, while 
the prediction of the electromagnetic output (and the stellar mass density to some extent) 
do depend on assumptions concerning the stellar IMF. The observational measurements are 
done for the electromagnetic output and SFR is derived by assuming an IMF, therefore the 
observational estimates of SFR suffers from the uncertainty in the assumed IMF. Further 
uncertainty is expected from dust extinction. This is particularly important when we deal 
with far UV hght. 

We also remark that the star formation rate in the simulations could depend on the 
details of the SN feedback prescriptions. If the feedback is turned off in the simulations, the 
heating of the gas is underestimated and therefore the global SFR density at low rcdshift is 
overestimated (i.e., ^/i^/szca/ overcooling of the gas). In the TVD simulation used in this paper 
and the SPH simulations on which the H&S model is based, the strength of the feedback was 
set based on the energetic argument of SN explosion and was also cahbrated against various 
observations such as the star formation rate of Lyman break galaxies, metallicities of the Lya 
forest and intra-cluster gas, and the galactic wind speeds observed in starburst galaxies. The 
current spatial resolution of the TVD simulation (physical cell size of 83/(1 + z)h~^ kpc) 
is not adequate to resolve the details of the inner structure of galaxies, but the overall 
star formation rate can be simulated with some confidence as the energetics are faithfully 
simulated on scales slightly larger than the real galaxies. 



4. StellEir mass density 

As we noted, the integral of SFR from t = to to is not equal to the stellar mass density, 
as the gas recycles into the interstellar medium. Figure 6a shows the evolution of the global 
stellar mass density as a function of redshift. The model results take account of the 
time-varying recycling fraction derived from BClibOS. All results assume the Chabrier IMF 
(0.01 - 100 M©). 

The values of at 2; = after the gas recycling correction are 0.0021, 0.0028, 0.0021, 
and 0.0023 for the Fossil, TVD, H&S, and SA models, respectively (see Table 3). The 
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empirical estimates converges to ^ (2.4 — 3.6) x 10~^ with the Chabrier IMF (0.01 — 
100 Mq), or (1.8 — 3.9) x 10~^ taking the errors into account. The model values are consistent 
with this range with a lower value favored. We remark that the conversion factor for from 
Salpeter IMF (mass range 0.1 - 100 Mq) to Chabrier IMF (0.01 - 100 Mq) is 1.4 for the 
same amount of luminosity output at M > IMq. The data arc taken from Cole et al. 
(2001, green filled circle), Panter et al. (2004, black open circle, the size of the claimed error 
bar is comparable to the symbol size), Fukugita & Peebles (2004, black inverted triangle), 
Kochanek et al. (2001, black open pentagon, slightly offset from z — ior clarity), and Bell 
et al. (2003, solid pentagon). 

The situation is different for the stellar mass for higher z. The figure compares the 
model curves with the observational analyses given by Rudnick et al. (2003, magenta filled 
squares. Table 2), Dickinson et al. (2003, cyan open stars) Brinchmann & Ellis (2000, red 
open squares), Fontana et al. (2004, blue open crosses. Table 4, 'observed' column), and 
Glazebrook et al. (2004, black filled triangles) all transformed to the Chabrier IMF. 

At z > 1 the observational estimates all fall significantly short of the model. Roughly 
speaking, all models predict that ~ 60% of the present stellar mass was formed hy z — 1, 
whereas the observation indicates it is as small as 20 — 30%, a gross disagreement, as also 
discussed by Nagamine et al. (2004) earlier. Wc emphasize that this is likely to be an 
observational problem, since the straightforward integration of the empirical star formation 
rate, corrected for the recycling factor, yields a stellar mass substantially larger than what 
is observed, indicating a gross underestimate of the stellar mass in high-rcdshift galaxies in 
currently available analyses. We note that possible errors in extinction corrections for z > 1 
galaxies do not solve this problem. 



5. Metal mass density 

Given the star formation histories, one can compute the metal mass density as a function 
of redshift as 

= Ypz, (7) 

where wc take Y = 0.023 using the prescription of Arnctt (1996, Section 14.4) adjusted to 
the Chabrier IMF. (This value corresponds to the Scarlc-Sargent yield Y/{1 — R) = 0.042.) 
The evolution of metal mass density is shown in Figure 7 together with empirical estimates. 

At 2; = 0, all models agree well with the estimate by Fukugita & Peebles (2004), Qz = 
(0.8 ± 0.25) X 10^"^ (eq.[94]), which excludes the contributions from stellar remnants but 
includes those from warm galactic halos. Dunne et al. (2003, Table 1) give an estimate 
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higher than Pukugita & Peebles (2004) by 50%, and the TVD model is consistent with their 
central value. Other model predictions are consistent with Dunne et al.'s estimate at its 
lower edge. 

Dunne et al. (2003) also give an estimate of the metal abundance at 2; = 2.5 including 
damped Lya systems. The models are consistent with their estimate at z — 2.5, although 
the TVD and H&S model lie somewhat above the given error bar. Bouche et al. (2006) 
and Prochaska et al. (2006) presented metal abundance from high-z galaxies and damped 
Lya systems lower than Dunne et al. They are, however, taken as lower limit, first because 
of missing contributions of faint galaxies as discussed by the authors, but perhaps more 
importantly by the omission of metals in warm galactic halos, which may contribute by 
additional 50-100%. Therefore we consider that the "missing metal problem" is not really a 
problem. 



The last test concerns the consistency with the energy in the radiation field, the product 
of stellar evolution. To estimate EBL we use the compilation of Hauser & Dwek (2001) for 
the source of observations. We write the integral of the flux over the range between Ai, and 
A2 (in //m), as 



in units of nW m"^ sr~^. Por the optical to the near IR region direct observations of 
extragalactic background yield /ebl[0.16, 3.5] ~ 60, whereas the integration of galaxy counts 
give /ebl[0.16, 3.5] ~ 18. We may take these two numbers as upper and lower limits, since 
the former is derived from difficult observations that are apt to be contaminated with local 
emission, while the latter would miss light from the outskirts of galaxies. Por the far IR that 
is dominated by dust we adopt Hauser-Dwek's estimate: /ebl[3.5, 140] = 11 — 58, where the 
upper value comes from fluctuation measurements and the lower from integration of resolved 
sources. Measurements are more secure for the submillimeter region: /eel [140, 1000] = 
15 ± 2. Adding the three together, we take the total obscured EBL flux to be 



6. Energy in the radiation 




(8) 



/eel [0.16, 1000] = 42 - 135 nWm'^sr 



-1 



(9) 



which ranges over a factor of 3. We avoid quoting the central value, since this range represents 
predominantly systematic uncertainties and what value should be taken as central is rather 
a matter of interpretation. 
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The bolometric EBL and the comoving bolometric luminosity density at redshift 
jhoi{z)-i s-re related by 



'EBL 



9.63 X IQ-^h-^ 



dt 



dz 

oo 







dz 
l + z 

Jbol(^) 

LqMpc" 



Ho 



dt 



dz 



dz 
l + z 



where 



dt 



dz 



[i + z)^nu{i + zf + nK 



(10) 



[nWm-^sr-^], (11) 



(12) 



for a flat A universe, and jboi(-2^) can be obtained by 

jboi(i)= / p*(r)[Lboi(t-r)/M]dr, (13) 

JO 

with p*(t) the comoving SFR density in units of MQyx~^ Mpc~^, and L\,o\{t)/M (erg Mq^) 
the bolometric luminosity per mass of a stellar population with age t. 



6.1. Bolometric luminosity density 

Figure 8a shows the comoving bolometric luminosity density as a function of redshift 
for the Fossil, TVD, H&S, and SA models. For the latter three models a solar metallicity is 
assumed. Figure 8b shows the spheroid and disk components separately for the Fossil model. 
We use the bolometric electromagnetic output to minimize uncertainties due to corrections 
for dust obscuration. 

We take Kashlinsky (2005) 's summary for observed luminosity density at the present 
epoch: jboi(0.2-2/xm) = (9.8±1.2) x 10^^ /lergs"^ Mpc~^ andiboi(12-100/xm) = (1.5±0.3)x 
10^^ /lergs"^ Mpc~^. Adding these two and using h — 0.7, we obtain jboi(0.2 — 100//m) = 
(3.98±0.43) x lO*L0^boi Mpc~^. Note that the energy density in the range of the gap, 2 — 12/im 
is expected to be small from the observations of the EBL. We expect the contribution from 
A > 100 /im as /ebl[120, 1000] = 10-15 nW m'^ sr^^ (Hauser & Dwek 2001; Kashhnsky 
2005). This yields jboi(0.2 - 1000/xm) = (4.4 ± 0.5) x lO^Le^boi Mpc"l Being conservative 
on the error, we adopt 

logjboi = 8.6±0.1, (14) 

as our estimate of the bolometric luminosity density at 2; = 0, which is shown in Figure 8. 
This estimate agrees with that of Bell et al. (2003, see their Fig. 15). It is interesting that 
there is no discrepancy amongst the models and between the models and observations at 
z — while the model predictions are rather divergent to one order of magnitude at high 
redshift. 
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6.2. 



EBL 



From the bolomctric luminosity density shown in Figure 8, we compute the total EBL 
using equation (11) as a function of redshift for each model of cosmic star formation history. 
Our calculation does not include the contribution from AGN, which amounts to 7%, as esti- 
mated in Nagamine et al. (2006, in preparation). The model predictions are all convergent 
on 7ebl ~ 37 — 51 nW m'^ sr~^ at 2; = 0, which are to be compared with the observation, 
-^ebl[0.16, 1000] = 42 — 135 nW m~^ sr~^ quoted above. This indicates that the models 
are consistent with the observation only at its lower edge, suggesting that the current EBL 
observations might still be contaminated from non-cosmological sources, whereas the inte- 
gration over resolved sources yields the value theoretically expected ^. The model prediction 
would even undershoot the observations if the Salpeter IMF were adopted. 

We may have another look at this problem by considering the ratio of the EBL to the 
bolometric luminosity density at ^ = 0. We define the dimensionless parameter 



where tn = ^/Hq is the Hubble time. This parameter compares the EBL with the luminos- 
ity density multiplied by the Hubble time, and the advantage is that it is independent of 
uncertainties in the stellar IMF provided that it does not vary as a function of time. From 
Figures 8 and 9, we compute 77^^, as shown in Figure 10. The results are indicated by arrows, 
and the solid line shows 77* as a function of r for the delayed exponential models. They are 
compared with the observation indicated by shades that ranges from 77^ = 0.61 to 3.11 from 
equations (9) and (14). All models result in 77^ = 1.0 ±0.2, which is consistent with only the 
lower part of the apparently observed range. Therefore, the problem is IMF independent. 



The evidence we have discussed in the previous sections indicate that we now have 
reasonable understanding of the cosmic star formation history. There are, however, two 
items that do not fit well to the scenario. They are the stellar mass density at 2; > 1 and 
the EBL from direct measurements of sky. 

We showed that the observationally derived stellar mass densities at 2; > 1 are signif- 
icantly smaller than the prediction of all the models, while the models agree well with the 



^See Note Added in the end of this paper. 
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7. Inconsistency among observations 
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observed SFR up to 2; ~ 5. This means that the observational data of SFR and stellar mass 
density are inconsistent with each other. The integration of SFR, model or observation, gives 
the correct stellar mass density at z = 0, which is a strong constraint. Hence, this is clearly 
an observational problem that observations underestimated the stars locked in galaxies at 
z > 1. Is it possible that the estimates of the stellar mass density are correct at all redshift 
and the SFR is mis-estimated? This is unlikely, because, to give the curve of the observed 
stellar mass density at 2; > 1, SFR must decline sharply from z 0.5 to higher redshift 
while the SFR density z < 0.5 is substantially higher than is measured. Therefore, we may 
ascribe the problem to underestimates of the stellar mass density at 2; > 1, presumably from 
poor understanding of the optical/NIR luminosity functions for sub-luminous galaxies.. 

The observation of EBL harbors another problem. All model results are consistent with 
the EBL only when we take the value obtained by integration of resolved sources, in both 
optical and infrared bands. In the optical band, direct measurements (Bernstein et al. 2002) 
give the EBL 3 times that from integration of resolved sources (Madau & Pozzetti 2000; 
Totani et al. 2001). This factor becomes 5 in the infrared (see Hauser & Dwek 2001, and 
references therein). The directly measured EBL in either of the cases gives too bright a 
flux to be accounted for. This conclusion does not depend on cosmological theories, nor on 
the IMF. It is a result of requiring consistencies between stellar mass density today and the 
light that has been emitted by those stars with subsidiary information concerning the star 
formation history inferred from the Madau diagram. Therefore this is also an observational 
problem possibly due to insufficient subtraction of the Zodiacal light (see also Kneiske et al. 
2002; Mattila 2003) ^. Finally there is an inconsistency between the observed value of the 
luminosity density and the observed EBL for any plausible star formation history. 

One may argue for the possibihty that EBL in the infrared is actually bright as a result 
of emission from high-redshift stars which are not included in our modeling. We consider 
that this is unlikely, since we do not find at 2; = extra stars that were responsible for 
such emission at high-redshift, and the amount of stars at 2; = is well documented. If 
there were very massive stars such as Population III (e.g. Cambrcsy ct al. 2001; Matsumoto 
2001; Cooray et al. 2004; Kashlinsky et al. 2004; Kashlinsky 2005; Matsumoto et al. 2005), 
these objects would die quickly and will not be counted as stars at the present epoch. How- 
ever there are several difficulties with this hypothesis, such as much higher SFR at 2; > 7 
than at 2; = 3 — 5 (Fernandez & Komatsu 2005), possibility of over-enriching the inter- 
galactic medium by metals (Salvaterra & Ferrara 2003; Ricotti & Ostriker 2004), physically 
unreahstic absorption-corrected spectra of distant TeV blazars (Dwek et al. 2005b, a), and 



See Note Added in the end of this paper. 
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overproduction of soft X-ray background (Madau & Silk 2005). 

8. Conclusions 

We have considered three physically motivated approaches: one with ab initio theoretical 
models that are constrained basically by the observation of WMAP a,t z ^ 1000, another with 
the Fossil model that uses theory of stellar evolution constrained by the SDSS observation 
at z ~ 0.1, and yet another purely from observations at z — — 6. We considered various 
physical quantities that witness the star formation history: stellar mass density, luminosity 
density, bolomctric luminosity density and EBL, and wc found a general agreement among 
results from these three approaches within the errors the current observations indicate. In 
particular, wc found that the exceedingly simple Fossil model having two populations with 
parameters adjusted to the SDSS observation — an update of the Scarle et al. (1973) model 
— provides a very effective description of the star formation history predicted by ab initio 
theoretical models, and it also exhibits a good fit to the suite of tests we have apphed. 

The only exception was in the stellar mass density at high redshift, where not only 
the Fossil model but all models predict higher values than the observationally derived ones. 
The observational results do not look consistent with what are inferred from the empirical 
Madau plot. This is not a matter of normalization of the star formation histories shown in 
Fig. 5, since the data on stellar mass density at z — are now well converged, and they are 
consistent with the integral of the SFR when the gas recycling is taken into account (see also 
Fukugita & Peebles 2004). If the dust extinction effect were significantly weaker than what 
we adopted here and the normalization of the SFR history were brought down in such a way 
that the integral of the SFR is consistent with the current data on stellar mass density at 
z > 1, such a SFR history would be inconsistent with the measured stellar mass density at 
z — 0. Therefore, we regard the disagreement in the stellar mass density at z > 1 between 
the models and observations as an observational problem. 

We also noted that the predictions of the EBL from the models all converged at the 
lower edge of the best current observations, giving Jebl = 43 ± 7nWm~^sr~^ This implies 
that the true EBL flux is presumably close to the value obtained by integrating over resolved 
sources. 

Before proceeding further, we should ask if an incorrect modeling of the stellar IMF 
could be the culprit. The number and mass in stars less massive than O.4M0 are quite un- 
certain as they do not significantly affect observed spectral properties and their contribution 
to dynamically determined mass estimates is obscured by the uncertainties in the dark mat- 
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ter component. For our analysis we consistently adopted the Chabrier IMF (0.01 — lOOM©). 

Let us now suppose that, due to a higher than expected fraction of low-mass stars (as com- 
pared to the Chabrier IMF), the effective mass-to-light ratio values were all to be increased 
by a factor of yU (yU ^ 1.6 for the Salpeter IMF; values of < 1 are also permitted if there 
were a even stronger turn-over of the stellar number density with decreasing mass as com- 
pared to the Chabrier IMF). This change would alter the theoretical models by reducing the 
EBL by a factor of about ji. The observed values of EBL would of course be unchanged. 
Thus, increasing the contribution of low-mass stars as compared to the Chabrier IMF would 
exacerbate the problem of under-predicting the EBL for all the models. On the other hand, 
reducing the contribution of low-mass stars even more so than the Chabrier IMF would ease 
the disparity between the predicted and observed EBL. But then the models would overpre- 
dict the present-day luminosity density in various bands and disrupt the rough agreement 
with the observed values that we found in Table 3. Therefore no simple adjustment of the 
stellar IMF can be made to remove all discrepancies. In addition, the parameter rji, was 
constructed (Equation [15]) in a fashion so as to avoid uncertainties in the stellar IMF and it 
too shows a discrepancy between observed and predicted value. Therefore it seems inevitable 
that the models of cosmic star formation history that we considered in this paper cannot 
account for a total optical-to-IR EBL of ~ 100 nW m~^ sr~^. We would suggest that the 
direct observation of the EBL in sky may still be contaminated by local emissions. If future 
observations would converge to a higher EBL value, as the current direct EBL observations 
indicate, we would get into a trouble in modeling the history of galaxies: we may have 
to invoke additional components, such as mini black holes or population III stars with a 
top-heavy IMF that are not included in the current modeling as discussed in Section 7. 

An interesting feature of our Fossil model is that it predicts that the spheroid component 
formed predominantly at z > 1.5, and the disk component formed mostly at z < 1. The 
early formation of the spheroidal component is consistent with some recent direct numerical 
simulations (e.g., Abadi et al. 2003; Naab et al. 2005). This agrees with the observational 
knowledge known for some time (e.g., Fukugita et al. 1996) and the early work by Searle et al. 
(1973) and predecessors. Although the current ab initio models cannot distinguish between 
these two components, the Fossil model provides one example how the total population can 
be separated into the two population, consistently with observations. As a consequence of 
the early formation of the spheroidal component, our models clearly indicate a relatively 
constant star formation rate from z = 7 to z = 2. It will be interesting to see if this 
prediction is confirmed or refuted by ongoing observational programs. 

We also note that the rapid decline of the spheroid formation from 2; = 3 to 2; = 1 in 
the Fossil model works in favor of having a population of Extremely Red Objects (EROs; 
e.g. McCarthy 2004) at ^ = 1 as found by recent observations (e.g. Cimatti et al. 2004; 
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Glazebrook et al. 2004; McCarthy et al. 2004). The existence of these EROs were first 
regarded as a challenge to the hierarchical CDM models (Somerville et al. 2004), but was 
later shown (Nagamine et al. 2005b, a) that the overall space density of EROs could be 
accounted for in cosmological hydro-simulations. Having a high value of SFR density at 
high-redshift {z > 4) as our models (except the SA model) helps to resolve the issues with 
the existence of high-redshift EROs. 

We conclude that the agreement among the three different approaches, as summarized 
in Table 3, is encouraging, and that the Fossil model provides a simple effective description of 
cosmological star formation history with very few free parameters. Furthermore, the general 
agreement between the TVD, H&S, and the Fossil model suggests that the cosmological hy- 
drodynamic simulation based on a cold dark matter model is providing a reasonably accurate 
picture of cosmological star formation history without a fine tuning of input parameters. 
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with the SFR data points in electronic form. This work was supported in part by grants 
NAG5-13381, AST 05-07521 and NNG05GK10G. MF received support by Grant in Aid of 
the Ministry of Education at the University of Tokyo and by the Monell Foundation at the 
Institute for Advanced Study. The cosmological TVD hydrodynamic simulation was carried 
out at the National Center for Supercomputing Applications (NCSA). 



NOTE ADDED: 

After the submission of our paper, Aharonian et al. (2006) published an important upper 
limit to the NIR background from the HESS observation of high energy gamma rays from two 
blazers, H2356— 309 and lESllOl— 232. The absence of absorption of gamma rays excludes 
the high value of EBL inferred from NIR fluctuation measurements, while the EBL from 
the source count is consistent with the observation. This agrees with our conclusion derived 
from the analysis given in Sect. 5.2. We also note that the source count recently measured 
with Spitzer (Dole et al. 2006) agrees with the EBL we concluded. 



-19- 



A. Mean metallicity of bulge 

We calculate the mass-weighted mean metallicity of the early type galaxies and the 
bulge component of spiral galaxies by 

, ^ ( Z(M)M %dM , ^ 

(z) = ^ — ; — Ai 

^ JZiL)L{f)^J[dL 

/L-^^dlogL ' ^ ^ 

where Z{M) and Z{L) arc the metallicity of the bulge component as a function of stellar 
mass M and luminosity L. In equation (A3) we used the scaling (M/L) oc L*^'^^ for the 
stellar mass-to-hght ratio of late- type galaxies (Bernardi et al. 2003; Vale & Ostriker 2004). 
We estimate the function Z{L) from the relation between metallicity and velocity dispersion, 

[Z/H] = 0.53 (log (7 - 2.173) + 0.15 (A4) 

obtained by Nelan et al. (2005). Using the Faber- Jackson relation (Faber & Jackson 1976), 

220 / L ^ 



^1 I km s~\ (A5) 
we obtain 

Z{L) = 0.13 log (^-^^ + 0.16. (A6) 

We assume the same relation for bulges of late type galaxies. 

We take r-band luminosity functions (LFs) of different morphological types of galaxies 
from Nakamura et al. (2003), and scale down the luminosity of SO/a-Sb, Sbc-Sd galaxies by 
factors of of 0.40 and 0.24 (see Table 1), respectively, to exclude the contribution from the 
disk component. The result of the integral (A3) is presented in Table 4. 

A caveat is that Equation (A6) may not be entirely appropriate for the bulges of late 
types galaxies, because the relation suggests the metallicity of the bulge of the Milky Way 
would be super-solar (0.5 dex), while observations suggest sub-solar metallicity (-0.5 dex) 
for bulge stars of the Milky Way (e.g., Freeman & Bland-Hawthorn 2002). We may have 
overestimated the mean metallicity for bulges of disk galaxies. 
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B. 



Determination of the Normalization of the Fossil Model 



The disk and the bulge components of the Fossil model are normalized to the observed 
luminosity densities jK,tot = 4.1 x 10^ hjoLQ^xMpc'^ (Bell et al. 2003; Cole et al. 2001) and 
jB,tot = 1-3 X 10^ /i7ol/0,s Mpc^^ (Fukugita & Peebles 2004). The luminosity densities are 
decomposed into spheroids and disk components, as 



where the subscripts b and d stand for the spheroid and the disk, with the coefficients the 
inverse of the stellar mass-to-light ratios bi — {Lk/M^)i„ 62 = {Lb/M^)i„ di = {Lk/M^)^, 
and d2 = {Lb/M^)^. The BClib03 model give 61 = (0.86)-^ = 1.16, 62 = (5.59)-^ = 0.18, 
di — (0.55)~^ = 1.82, and ^2 = (1.20)~^ = 0.83. We solve the equations and obtain 
p^,;, = 1.63 X 10^ M© Mpc"^ and p^,^ = 1.21 x 10^ M© Mpc"^. The corresponding stellar 
mass density parameters are = 0.00120 (bulge) and 0.00089 (disk) as presented in Table 2 
and 3, which means the mean bulge-to-disk mass ratio (Mh/Md) = 1.35 (or equivalently 
Mb : Md = 57% : 43%). This bulge mass fraction is somewhat smaller than the one obtained 
by explicit summation over the galaxy sample, but it is within the expected uncertainty of 
~ 20%. We note that this method does not work well if one uses the r-band luminosity 
instead of S-band luminosity, because r-band is not sensitive enough to young stars that 
dominate the disk component. 
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Table 1. Bulge-disk luminosity decomposition based on SDSS data 





E 


SO 


SO/a-Sb 


Sbc-Sd 


Irr 


jr 


Jr 


0.21 


0.41 


1.00 


0.37 


0.02 


2.00 


B/T^ 


1.00 


0.62 


0.40 


0.24 


0.00 


0.95 


D/T^ 


0.00 


0.38 


0.60 


0.76 


1.00 


1.05 



^ r band luminosity density contributed by each Hubble type of galax- 
ies from Nakamura et al. (2003). 

^''^ Bulge-to-total and disk-to-total luminosity ratios for each type of 
galaxies from Ohama (2003). 

Note. — The last column gives the r band luminosity densities 
for the bulge and disc components in units of 10^ /iLq^^ Mpc~^, which 
yield jr{B) / jr{D) — 0.91. Extinction corrections are applied to the 
disk component, using Ar — 7ijlog(a/6), where a/b is the axis ratio, 
and 7r = 1.15 + 1.88(log ^^20% - 2.5) = 1.24 at log 1^20% = 2.55 for 
R* = -21.09 from Sakai et al. (2000, see also Tully et al. (1998)): for 
log (a/b) = 0.22 we obtain (Aji) = 0.27, where the difference between R 
and r is ignored. This leads us to the corrected bulgc-to-disk ratio for 
the luminosity density, >(5)/>(L))|corr = >(5)/>(L')10°-^<^«) = 0.71. 
We assign a 20% error to this quantity. 
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Table 2. Properties of the Fossil model at 2; = 



Population 








u-g 


g-r 






Bulge 


0.00120 


1.5 


1.5 10.5 


1.78 


0.84 


5.59 3.19 


0.86 


Disk 


0.00089 


4.5 


0.8 7.0 


0.94 


0.43 


1.20 1.18 


0.55 



^ Characteristic time-scale in units of Gyr [defined in Equation (1)]. 
^ Metallicity {Zq = 0.02). 

^ Mean age of stars in units of Gyr. 

Stellar-mass-to-light ratio in solar units for the i?, r, X-bands and the Chabrier IMF 
(0.01 — 100 M0). Stellar masses include the remnants. 
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Table 3. Physical quantities at z = for the models used in the text 



Model 








Iebl^ 




iboi'' 






Jr 




u- g 


g -r 


Fossil 


0.0021 


2.26 


0.58 


42 


0.86 


3.0 


1.4 


1.3 


1.5 


4.1 


1.54 


0.67 


bulge 


0.00120 


1.28 


0.36 


22 


0.022 


0.9 


0.2 


0.3 


0.5 


1.9 


1.78 


0.84 


disk 


0.00089 


0.99 


0.22 


20 


0.837 


2.1 


1.2 


1.0 


1.0 


2.2 


0.94 


0.43 


TVD 


0.0028 


3.06 


0.80 


51 


1.02 


3.2 


1.5 


1.4 


1.8 


5.1 


1.20 


0.61 


H&S 


0.0021 


2.23 


0.56 


37 


1.30 


3.2 


1.7 


1.5 


1.6 


4.1 


1.00 


0.50 


SA 


0.0023 


2.56 


0.61 


47 


1.37 


3.7 


2.0 


1.7 


2.0 


4.9 


1.04 


0.52 


Consensus 


0.0023 


2.52 


0.65 


43 


1.06 


3.1 


1.5 


1.4 


1.6 


4.4 


1.25 


0.59 


Model^ 


±0.0004 


±0.47 


±0.13 


±7 


±0.22 


±0.2 


±0.2 


±0.1 


±0.2 


±0.6 


±0.27 


±0.09 


Observed 


0.0018 






42- 


0.5 


3.2 


1.4 


1.2 


1.27 


3.4 






range 


-0.0039 






135 


-1.6 


-5.0 


-1.8 


-2.0 


-1.33 


-5.5 







^ Stellar mass density including the remnants for Chabricr IMF (0.01 — 100 Mq). 

^ Stellar mass density without the remnants (in 10* Mpc~^). 

^ Stellar remnant mass density (in 10* Mq Mpc~^). 

^ EBL intensity (in nW m~^ sr""-*^). 

§ SFR density (in lO-^Moyr"^ Mpc"'^). 

^ Bolometric luminosity density (in 10* boi Mpc~^). 

Luminosity density in U, B,r, K -hands (in 10* Lq^b Mpc~^). 
^ Average of the Fossil, TVD, and H&S models, with errors being the dispersion among the three models. 

Note. — Bolometric magnitudes Mb,q = 5.48, Mk,q = 3.28 and M^,© = 4.76 are adopted for the Sun. 
Model luminosity densities are compared with observations: js = (1.7 ± 0.3) x 10* Lq^b Mpc^^ (Nagamine 
et al. 2001a), Jb = (1.3 ±0.1) x 10* L0,ij Mpc"^ (Fukugita & Peebles 2004), > = (1.3 ±0.03) x 10* Lo^^Mpc"^ 
(at z = 0.1, Blanton et al. 2003a), Jk = 'i-ltH x 10^ Lq^k Mpc'^ (BeU et al. 2003), (5.0±0.5) x 10* Lq^k Mpc'^ 
(Kochanek et al. 2001), jxs = (4.0 ± 0.6) x 10* Lq^ks Mpc"^ (Cole et al. 2001). All values assume h = 0.7. 
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Table 4. Mean metallicity of spheroids 



lype {Z/Zq) - jr [10« h Lr,e Mpc"^] ^ 



E-SO 1.44 0.70 

SO/a-Sb 1.15 0.39 

Sbc-Sd 1.15 0.082 

ibtai Tss rrr 



Mass weighted mean metalhcity. 
^ r-band luminosity density from the spheroid component. 




Fig. 1. — Star formation histories of delayed exponential decay model as functions of cosmic 
time {panel a) and redsliift {panel b). The normalization is fixed to equation (2) for each 
model. The lines are for r = 0.1 Gyr (black long-dashed), 1 Gyr (magenta dotted), 2 Gyr 
(blue solid), 3 Gyr (cyan dashed), 5 Gyr (green dot-dashed), and 10 Gyrs (red long-short- 
dashed). The blue thick solid line is the composite two population model with r = 1.5 Gyr 
and 4.5 Gyr, which is taken as our 'Fossil' model. 



- 32 - 




Fig. 2. — Gas recycling fraction i? as a function of cosmic time (and redshift), computed 
using BClib03 for solar metallicity and the Cliabrier IMF (0.01 - 100 Mq). 'SSP' stands for 
'Simple Stellar Population' with an instantaneous burst at t = 0. 
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Fig. 3. — Panel (a): Color-color diagram 'm u — g vs. g — r plane. The three lines represent the delayed 
exponential decay models with Chabrier IMF for metallicities: Z/Zq — 0.2 {black open squares), 1.0 (blue 
solid triangles), and 2.5 {red open triangles), with t running along the line from top to bottom as indicated 
in the legend. The contour shows the color distribution for SDSS galaxies taken from Blanton et al. (2003b). 
Panel (b): Same as (a), but with the histogram for the color distribution of SDSS galaxies projected onto 
the axes together with Gaussian fits. The two boxes indicate the location of the peaks of the histogram, 
with the widths twice the dispersions of the Gaussian. The asterisks indicate the bulge component (with 
r = 1.5 Gyr, Z/Zq = 1.5) and the disk component (with r — 4.5 Gyr, Z/Zq — 0.8) of the 'Fossil' model. 
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14 




5 10 



T [Gyr] 

Fig. 4. — Mean age of stars as a function of decay time r for the delayed exponential model. 
The dotted lines shows the values of r (r = 1.5 Gyr and 4.5 Gyr) adopted for the 'Fossil' 
model. The arrows indicate the shift in the mean age when stellar remnants are removed 
from the calculation, while the curve includes the contribution from the dead stars that are 
not emitting light at the present time, therefore push the age up. 
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Fig. 5. — Star formation rate density as a function of redshift. The curves represent the model 
predictions specified in the legend. The data are taken from (from low to high redshift): 
Heavens et al. (2004, 3 asterisks at z ~ 0), Nakamura et al. (2004, open inverted triangle at 
z = 0), Lilly et al. (1996, open circles), Norman et al. (2004, filled triangles), Cowie et al. 
(1999, open diamonds), Gabasch et al. (2004, open squares), Reddy et al. (2005, cross at 
z = 2), Barger et al. (2000, open pentagons at z = 2 and 4.5), Steidel et al. (1999, filled 
pentagons at z = 3,4), Ouchi et al. (2004a, filled circles at z = 4,5), Giavalisco et al. (2004, 
open triangles at 2; = 3 — 6), Bouwens et al. (2005, filled square at z = 6), and Thompson 
et al. (2006, open stars without error bars). The data are converted to the values with the 
Chabrier IMF and common values are assumed for dust extinction for the UV data. See text 
for details. 
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Fig. 6. — Panel (a): Stellar mass density (including remnants) as a function of redshift and 
age. The top axis indicates the age of the Universe. All data assume the Chabrier IMF 
and h = 0.7. The source of data are given in the text. Panel (b): Growth of the stellar 
mass density normalized by the value at z = 0. For the Fossil model, the bulge and disk 
components are also shown separately. 



-37- 



"1 r 



"T r 



Fossil 

TVD N1024L85 
H&cS 



I — I 

CO 



O 




o 







6 



redshift 



Fig. 7. — Metal mass density as a function of redshift. The curves represent the model 
predictions specified in the legend. The data are taken from (from low to high redshift): 
Dunne et al. (2003, red bars at 2; = & 2.5), Fukugita & Peebles (2004, blue cross at z = 0, 
shifted for clarity), Bouche et al. (2006, magenta lower limit ai z = 2, for galaxies), and 
Prochaska et al. (2006, black lower limit at z = 2, shifted for clarity. Only for damped Lya 
systems and super Lyman Limit systems). 
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Fig. 8. — Panel (a): Comoving bolometric luminosity density as a function of redshift for 
the models indicated in the legend. Panel (b ): Same as (a) for the Fossil model, decomposed 
into the two components. The data at ^ = (slightly offset to a positive z value for clarity) 
is the observational estimate given in Equation (14) 
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Fig. 9. — EBL as a function of redshift for the models shown in the legend. The observa- 
tionally allowed range at z = 0, equation (9), is indicated by the error bar. 
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Fig. 10. — Parameter rj^, of Equation (15) computed for the models used in this paper (shown 
by arrows). The variation of rj^, is also shown as a function of decay time-scale r for the 
delayed exponential model of star formation. The shaded region shows the observationally 
allowed range. 



